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Abstract 

Starting from the diffusion equation at beta random matrix hard edge obtained by 
Ramirez and Rider (2008), we study the question of its relation with Lax pairs for 
Painleve III. The results are in many respects similar to the ones found for soft edge 
by Bloemendal and Virag (2010). In particular, the values beta = 2 and 4 (but not 
beta = 1) allow for a simple connection with Painleve III solutions and Lax pairs. 
However, there is an additional surprise for a special relation of parameters where a 
simple solution of the diffusion equation can be obtained, which is a one-parameter 
generalization of Gumbel distribution. Our considerations can be extended to the other 
Painleve equations since the corresponding diffusions are in fact known as nonstationary 
(imaginary time) Schrodinger equations for quantum Painleve Hamiltonians. We also 
track the hard-to-soft edge limit transition in terms of our Lax pairs. 
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1 Introduction 



Beta-ensembles of random matrices were originally defined by their eigenvalue probability 
distribution depending on real nonnegative parameter denoted /5 by Dyson [15j who first 
introduced the probability density function (p.d.f.), 

1 

Pf^ih, ■■■,ln) = ^-l[\h-ljf- l[e-''^^'^\ (1.1) 

i<j k=l 

where V{x) is called the potential and is a normalization constant. This expression is 
the joint eigenvalue p.d.f. for three important invariant ensembles of random matrices, whose 
entries have probability distribution symmetric with respect to the real orthogonal group 
(/9 = 1), complex unitary group (/3 = 2) and symplectic group (/3 = 4). For the other values 
of /3 there had been no known random matrix model realizing the above eigenvalue distribution 
until 2002. Then Dumitriu and Edelman [H] discovered tridiagonal random matrix ensembles 
which lead to (1.1) for V{x) = (Gaussian /3 ensemble) and V{x) = —x + const. ■ \nx 
(Laguerre (or Wishart) /3 ensemble). This discovery spurred a lot of further research on /3- 
ensembles and their large size n limits. New and different matrix models for /3-ensembles, 
e.g. invariant rather than tridiagonal [2], have been also constructed. Even without the 
corresponding matrix ensembles, the eigenvalue distributions (1.1) have physical significance 
as description of Coulomb gas(fluid) restricted to a line, which has application in problems of 
conductivity and other transport phenomena in various disordered systems, see e.g. [19] and 
references therein. 

Their theoretical importance is also due to the relation with conformal field theory exposed 
by Virasoro subalgebra constraints they satisfy [31 [D EH]- This relation implies that for 
general /3 these eigenvalue distributions and matrix ensembles are natural subject of quantum 
integrable systems theory [7l|TTl|T9]. Indeed, recently their direct connection with (imaginary 
time) canonically quantized Painleve Hamiltonians was shown for some special potentials [25] . 
There certain Fokker-Planck (FP) (or diffusion-drift) equations reappeared which had been 
found before by completely different, probabilistic, approach in [S] based on earlier work [2B] 
for random /3-matrix soft edge and in [261 EZ] for the hard edge, the two important scaling 
limits for large n. 

In the last references, the tridiagonal matrices of [14] and their spikecQ modifications 
were the starting point to rigorously derive random Schrodinger equations in the large n 
limit, conjectured earlier in [TB]. Then they were reduced to first-order stochastic ODEs (of 
Langevin type though nonlinear) by the Riccati transformation first used likely by [22j. Then 
the corresponding probability distributions satisfied the FP equations featuring the quantum 
Painleve II (for the soft edge) and quantum Painleve III (for the hard edge) Hamiltonians. 
In [8] the soft edge FP equation 

""^In statistics, a Wishart sample covariance matrix is called spiked if the underlying population covariance 
matrix has a single eigenvalue different from unity 
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(dt + ^dxx + it- x^)dx)J'l = 0, (1.2) 

where t is the scaled edge of the eigenvalue spectrum and x is the scaled spike parameter, 
and its bounded solution J^^ describing the largest eigenvalue probability for critically spiked 
large n Wishart /3-ensemble was shown to be related to Baik-Rains Lax pair ^ for (classical) 
Painleve II equation, however, only for /3 = 2 and 4. This pair. 



dx 



g ) \ -qx + q' x"^ - t - ) \ g 

"Mi'), (1.3) 



g \ q -X \ g 



is a modification of Flaschka- Newell [17j Lax pair found earlier to be relevant for the critical 
regime in the Tracy- Widom-to- Gaussian phase transition for the largest eigenvalue distribu- 
tion of n — 7- oo complex Wishart (/3 = 2) ensemble It was found in [S] that 

(t,x) = F2(t)/(t,x), Fiit.x) = Fl'\t){ip{t)f{t,x)+^{t)g{t,x)) (1.4) 

where -p2(t) is the Tracy- Widom distribution |3I], and, thus, for these special [3 only, the 
solution of the FP equation is proportional to either a component of the Lax pair eigenvector 
or a linear combination of its components. This turned out not so for the other values of /3 
and e.g., surprisingly, for [3 = 1, where results could have been expected similar to /3 = 4 
case, which is related to it by duality /3 ^ 4/(3 playing a big role in many aspects of /3-matrix 
models, see e.g. [531 [131 [21]- The part of the operator on the left-hand side of eq. (1.2) with x- 
derivatives is nothing but the canonically quantized Painleve II Hamiltonian which appeared 
as such in [25] (up to the change of sign of t). 

Remark. Just recently a nice formula for the spiked general /3 finite n Wishart matrix 
p.d.f. was obtained by Forrester [20] based on several earlier results for special cases: 

p,,aAh, . . . , /n) = ■ n - ^/ ■ n i^r'^'^'-'e-^^^/' / e^* n ( - 

where 6 is the (non-scaled) spike parameter and in statistical context a is the integer difference 
of dimensions of Gaussian rectangular nxm sample matrix X whose covariance matrix XX"^ 
is then nxn Wishart matrix. However, in general, the only restriction on a is a > —1 needed 
for convergence of the integrals over the above density. The hard edge limit arises here when 
n — )■ oo, n5 — )■ X, and the scaled spike parameter x is fixed. 

Another FP equation describes the hard edge large n limit of /3-ensembles, e.g. of (spiked) 
Wishart (Laguerre) ensemble, for the smallest eigenvalues restricted by hard wall of positiv- 
ity constraint. This was derived in [261 EZ] for different boundary conditions depending on 
whether the spiked or non-spiked case is studied. The equation is 
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(^d^ + ^x^d,, + (^{a + |)x - - e-^^ 9,^ ^^(x, fx) = 0, (1.5) 
and for the spiked case the boundary conditions are 

J^^(-oo,x) = 0, J^^(oo,x) = l, J^^(/i,0)=0. (1.6) 

The operator on the left-hand side of eq. (1.5) is this time 9^ plus the canonically quantized 
Painleve III (PHI) Hamiltonian which appeared as such in [25], up to the change of variable 
—7- —t. The FP equation (1.5) is the main object of study here, with the goal of ex- 
ploring the possible connections with classical integrability and comparing the results with 
that of soft edge Motivated also by [25j, we look for Lax pairs which give Painleve III 
equation as consistency condition such that the function J-"^ is a component of their eigen- 
vector. The dependence of such Lax pairs on spectral parameter, whose role is played by the 
spiking parameter x here, is rational and known in general for all Painleve equations |23] . 
With this example of PHI, we develop a general method of finding such quantum-classical 
correspondences, which should work equally well e.g. for the other Painleve equations and 
more broadly, a conceptual discussion of this approach to appear [30j. Here, essentially, the 
compatibility of FP equation (1.5) and a Lax pair for classical PHI, which brings additional 
algebraic constraints imposed on the elements of the Lax matrices, which not always can be 
successfully resolved to find nontrivial solution. When they can, some of them turn out to be 
first integrals for the ODEs involved which in effect prevents the whole system of ODEs and 
algebraic constraints from being overdetermined. 

One should mention also the studies of multi-spiked Wishart ensembles and their soft 
edge O [9] and hard edge [2^ where multidimensional generalizations of the above FP equa- 
tions emerge. Here again connections with quantum integrability are implied for general /3, 
see e.g. [2S] and references therein, and explicit determinantal formula for /3 = 2 distribution 
is known [1] in terms of the function / from eq. (1.3). But here we consider only the one-spike 
case. 



2 Main results 

We set up and solve the problem of finding possible classical Lax pairs for Painleve HI equa- 
tion, such that their first component satisfies eq. (1.5). The results obtained are in many 
respects similar to the ones for soft edge case |8j. Again, two values /3 = 2,4, but not (3 = 1, 
turn out to be the special ones under which there are the sought Lax pair related to Painleve 
HI. Besides these two, however, here appears another special value, more precisely, special 
relation: 

f + | = 2-a. (2.1) 

which admits a special unusual Lax pair depending on two arbitrary functions of fi. Actually 
we formulate our results in different variables. We change the independent variable fi for 
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t — e^, multiply eq. (1.5) by k, — /3/2, rescale variables as — )• t — >■ k,H and redenote 
a + 1/k — >■ a, and arrive at 

{Ktdt + x^d^^ + {ax -x^ - l/t)d^)J^^{t, x) = 0, (2.2) 
with the corresponding boundary conditions derived from eq. (1.6), 

jr^(0, x) = 0, ^"^(+00, x) = 1, J-^(t, 0) = 0. (2.3) 
Now we can write out the main outcomes explicitly. Consider Lax pairs of the form: 



d,(''")=B( ^"), (2.5) 



Then we obtain the following: 

Theorem 1. If j3 — 2, the solution J-'^{t,x) of eq. (2.2) is the eigenvector first component of 
the Lax pair of the form (2.4), (2-5) with 

<f> I ty4> \ 
tx'^ X I 

q^h 1 9^(9^-1) 1 kz^j^hzaj^l /' 



where function y{t) = /(p is a solution of the following Painleve III equation ('prime ' here 
and further on means derivative w.r.t. t), 

„ iy'Y 1 , q a , (a — 1) 1 „ s 

functions q{t) and y{t) satisfy the system of ODEs, 

ty = 2?^ - 1 - tS'^ - (a + l)t^y, (2.8a) 

and function h{t) is 

- 1 

h = —, (a-1). (2.10a) 

t^y ^ ' ^ ' 

If (3 — 2, the solution T^{t,x) of eq. (2.2) is also the eigenvector first component of the 
Lax pair of the form (2.4), (2.5) with 
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t2 t^y 



t^y4> f^x t 



where function y{t) = (j)'{t)/(f) is a solution of the following Painleve III equation 

functions q{t) and y{t) satisfy the system of ODEs, 

tV = 2g2 _ 1 _ ^3^2 + _ 2)t2^, (2.86) 



and i/iis time function h{t) is 



q2 - 1 

h = —^ — + a. (2.106) 

t y 



Remark. Function q from the Lax pair (2.6b) itself satisfies eq. (2.12) below, a variant 
of Tracy-Widom equation for /3 = 2 hard edge from [32]. If we denote it q{t] a) and denote 
7/(t; a) the solution of eq. (2.7b), then functions q and y from the Lax pair (2.6a) are nothing 
but q{t; 1 — a) and y{t; 1 — a), respectively. Nontrivial consistency conditions between the 
pairs (2.6a) and (2.6b) arise from the fact that they involve the same function in the 
eigenvectors. They lead to the following two algebraic identities, which can be directly verified 
from ODEs satisfied by q and y: 

yit- a)y{t- 1 - a) = + , _ = 0- 



The function q(t; a) also relates the solutions at /3 = 2 and /3 = 4. One has 

Theorem 2. If P = 4, the solution J^^{t,x) of eq. (2.2) is the eigenvector first component of 
the Lax pair of the form (2.4), (2.5) with 

l+q _a I 1— g S+ I 2tq'+a S+ . S+ 



2tx^ 2x ' 2 2ix2 ' 2(g2_l) x ' 2 

5- 2tq'~a I 1—q _a i_ l+q 

2tx^ 2(g2_l) !• "T 2 2te2 2x ~^ 2 



l+q S+ _|_ (2tg'+a)S+ 

i? = /.o(t)-/+| 2.--+^_y-i) 1, (2.11) 

2t^x 2t(g2-l) 2i2x 
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where function q{t) satisfies the following Painleve Ill-related equation, a variant of which 
first appeared in (but there (5 = 2): 

t{q' - = qitq'r + \q'iq' - 2) + Q - ^) g, (2.12) 

functions S+(t) and S^(t) are determined by 

S'+ _ qq' aq S'_ _ qq' aq q o _ i 2 (c. 

S+"g2_i + 2t(g2_i)' 2%2-i)' ^+^- " ^ ^' l^-^^) 

and function h{t) is 

1 f {2tqr - + ^ + a{a-2) \ _ ^^.u) 



" 2t \ 4(g2 _ 1) ' t ' 4 

Theorem 3. // eg. (2.1) holds, the solution J^^{t,x) of eq. (2.2) is the eigenvector first 
component of the Lax pair of the form (2.4), (2.5) with 

l+r 4> 
r — 1 1— r I ; 



l+r _|_ l+r I 

'4t2</)x 4Kt2<;i 24, 2t'^x 2Ki2 

For any functions r{t) and the Lax equations are satisfied identically with this pair. 

The FP eq. (2.2) is satisfied by the first component of the eigenvector only if eq. (2.1) is 
imposed. These conditions lead, however, to a simple explicit distribution function J-'^{t,x), 
unrelated to Painleve III, 

^^(t, x) = exp • (2.16) 

Thus, here only the partner of J-"^ , G{t,x), depends on the two arbitrary functions. It is given 
by 

G(t,x) = i^^^. (2.17) 

Only the above options for (3 give the possibility for J-'^(t,x) of eq. (2.2) to be related with 
classical Lax pair in this simple way. 

It is well known that when a becomes large, the smallest eigenvalues of random matrices 
separate from the hard wall at zero and the hard edge of the spectrum turns into a soft 
edge [in], see also [ISllISlEn]- Then, under proper scaling, 

X ^ a{l + a'^^^x), t ^ a'^{l + a'^/^t), « = |) (2-18) 
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the solution of eqs. (2.2), (2.3) transforms [2S1I2Z] into the bounded solution of eq. (1.2) 
after changing i^t — )■ t and kx x there. It is interesting also to follow this transition in 
terms of the above Lax pairs (2.6a), (2.6b) and (2.11). Doing this, we obtain the next 

Theorem 4. Under the scaling (2.18), when a — )■ oo, the Lax pairs (2.6a) and (2.6b) trans- 
form into the pair 

( ^'\_( -{qx + q') \( 

"^"^ )~\-qx + q' x^-t~q^ )\ 

where q{t) is the Hastings-McLeod solution of Painleve II equation, 

q" = tq + 2g^ (2.20) 

and function u{t) satisfies 

u = {q'f - tq^ - q\ u' = -q\ (2.21) 
The components of the eigenvector {J-'^,G) of (2.6a) scale as 

jrH ^ «-2/3j-^, G ^ e"'''*G^ (2.22a) 

while that of (2.6b) scale as 

jrH ^ a-^/'j^s, G ^ e°"''*G^ (2.226) 

Pair (2.19) is almost the same as (1.3), it becomes exactly (1.3) after changing the eigenvector 
as {J^^,G^) — 7- F2{t){f,g). Then u goes away from matrix B in (2.19) since u is also the 
logarithmic derivative of Tracy-Widom distribution F2(t): F2(t) = UF2. 

Theorem 5. Under the scaling (2.18), when a — )• 00, the Lax pair (2.11) transforms into the 
pair 




B^a'^'i ), (2.23) 

2 2 

where q is the same Hastings-McLeod solution of eq. (2.20), the functions S^,S_ satisfy the 
relations 

S'_^ = -qS+, S+S^ = 1, (2.24) 
and u is the same as in (2.21). The eigenvector components of (2.11) both rescale as 
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(J-^G) (2.25) 

Here our function S^{t) is seen to be the function denoted i?~^(x) in [8] (our t here is their 
x); in our notation, ip{t) = [S]/^ + S]/^)/2, ^{t) = [S]/^ - S]/^)/2, see eq. (1.4). The 
representation (2.23) for 13 = 4 case is still different from the representation in [S], here J-'^ 
and are linear combinations of / and g from (1.3). 

The rest of the paper is mostly dedicated to the derivation of the above results. In section 
3 we obtain the general set of ODEs and constraints relating the solution of FP equation 
(2.2) with a Lax pair of the form (2.4), (2.5). Then we consider different possible ways of 
constraint resolution in section 4. In section 5 we derive Painleve III (PHI) eqs. (2.7a), 
(2.7b) and Tracy-Widom eq. (2.12) from our equations. We then obtain the Lax pairs (2.6a), 
(2.6b) and (2.11) as well as Theorem 3 in section 6. In section 7 we consider the hard-to-soft 
edge limit transition and obtain Theorems 4 and 5. Some conclusions are presented in the 
last section 8. 

Acknowledgments. I would like to thank B. Rider for multiple useful discussions. 
The financial support by NSF grant DMS-0645756 is gratefully acknowledged. Also I would 
like to thank MSRI at Berkeley, CA, for support during my stay there in Fall 2010 there as a 
postdoctoral fellow in Special Program on Random Matrices and Random Processes, where 
my interest in the current problem was initiated. 

3 FP equation for the hard edge and Lax pair for PIII 

A Lax pair for Painleve III equation can be always found in the form such that 



«')+M)+i„(,). (3,) 



X 



and 



B = B,it)x + Boit) + (3.2) 

X 

as follows from singularity theory for ordinary differential equations (ODEs), see |23]. As 
follows from eq. (2.4), 

5... ( '^n)= (9.1 + L')( ^^), (3.3) 



G J ' ^ ' ^ V G 

and, using eq. (3.1), one finds that 



From eq. (3.1) one also finds 
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(ax - - l/t)L = + ""^^ ^^^^ + aLi -L^- LJt + iaL^ - LAx - Lqx\ (3.5) 

X'^ X 

Therefore the action of the hard edge FP operator from eq. (2.2) on the eigenvector of sought 
pair (2.4), (2.5) gives 

{Ktdt+x^d^^+{ax-x'^-l/t)d^) ^ = {KtB+x^{d^L+L'^) + {ax-x'^ -l/t)L) (^-^^ ^ ^ 



+ >.«a-. + {L..L.} + («-2)L,-W« ^ ^^^^ ^ ^ ^ 

X 

+(a - l)Li - L2 - y + (KiSi + {Li, Lo} + aLo - Li)x + (L^ - L^)x'^ ( "g^ ) ' ^^"^^ 

Let us denote the matrix on the r.h.s. of eq. (3.6) as T>h{L, B), then, for to be a solution 
of eq. (2.2), we must have 

where, according to the powers of x present in T>h{L, B), 

2 2 
ad{x,t) = J2 dk{t)x\ aa{x,t) = (3-8) 

k=-2 k=-2 

and the functions dk{t) and ak{t) are arbitrary so far. 

Besides the FP equation (3.6), which will give us some constraints, we have the Lax (or, 
more precisely, the zero curvature) equation: 

dtL^d,B + [B,Ll (3.9) 

which will give us the set of ODEs on the components of L and some more constraints. Since, 
by (3.1) and (3.2), 

[B_,,L2\ ^ [B_,,L,] + [Bo,L,] [B_,,Lo] + [Bo,L,] + [B,,L2] , 

-D, L\ = \ \ h 

X"^ X^ X 

+ [Bo, Lo] + [Bi, Li] + [Bi, Lo]x, (3.10) 
one obtains the following five matrix equations for the five powers of x present in (3.9): 
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x-'■.[B_^,L2]^0, (3.11) 

x-^ : dtL2 = + Li] + [Bo, L2], (3.12) 

: dtL, = Lo] + [So, Li] + [Si, L2], (3.13) 

x° : 9tLo = Si + [So, Lo] + [Si, L,], (3.14) 

: [Si,Lo] = 0. (3.15) 

Let us separate the trace parts of matrices L and S. i.e. write L — Lt ■ I + L, B — Bf I + B, 
where L and S are traceless. Then, eqs. (3.11) and (3.15) are equivalent to, respectively, 

S_i =S_u-/ + 6_(t)L2, (3.16) 

Si =Sif/ + 6+(t)Lo, (3.17) 

where 6+(t) and b_{t) are still arbitrary functions. The trace parts of the other matrix Lax 
equations give 

dtL2t = -S_u, (3.18) 

Lit — Cit — const., (3.19) 

dtLot = Si,. (3.20) 
Thus, eqs. (3.16), (3.18) and (3.17), (3.20) lead, respectively, to the expressions 

B_i ^ -dtL2t ■ I + b.{t)L2, (3.21) 

Bi = dtLot ■ I + b+{t)Lo. (3.22) 
The rest becomes the system of three matrix ODEs for traceless matrices: 

dtL2 = -6-L2 + [L2, 6-Li - So], (3.23) 

dtL, = (6_ - 6+)[L2,Lo] + [So,Li], (3.24) 

dtLo = b+Lo + [So - 6+Li, Lo]. (3.25) 
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Now we return to the consequences of FP equation (3.6), using also eqs. (3.21) and (3.22). 
One has five matrix equations for five powers of x present in T>h{L, B): 

x-^ : Ll - Ujt = d-2{I - Ts) + a-2<J-, (3.26) 
x-^ : Kt{b-L2 - dtL2t ■ I) + {^2, U} + (a - 2)L2 - Li/t = d_i(J - as) + a_i(j_, (3.27) 

x° : KtBo + Ll + {L2, Lo} + (a - l)Li - L2 - Lo/t = do{I - (73) + ao(7_, (3.28) 

x^ : Kt{dtLot ■ I + b+{t)Lo) + {Li, Lq} + aLo - = di{I - a^) + aia_, (3.29) 

x"^ : Ll - Lq = 62(1 - (T3) + 020-^. (3.30) 

The outcome of the last equations is 10 scalar relations (the rest just determine the 10 arbitrary 
functions on the r.h.s.). Among them there are 6 algebraic constraints on the components of 
L-matrices: 



{L2t + L2df - {L2t + L2d)/t + Ll^ = 0, 


(3.31) 


(2L2t - l/t)L2+ = 0, 


(3.32) 


{Lqi + Lq^Y — {Lot + Lod) + L1^ = 0, 


(3.33) 


(2Lot - l)Lo+ = 0, 


(3.34) 


6_ + 2Cu + a - 2)L2+ + (2L2t - l/t)L^+ = 0, 


(3.35) 


{Kth+ + 2Cu + a)Lo+ + (2Lot - l)L^+ = 0, 


(3.36) 



where we decomposed the traceless matrices as L = L^cr^ + L^ — L^crs + L^a+ + L_a^. The 
other equations are 2 ODEs for the traces: 

KtdtL2t = (Ktb- + 2Cit + a - 2)L2d + (a - 2)L2t + {Cu + Lu){2L2t - l/t) + {L2, Li}, (3.37) 



KtdtLot = -{Kth+ + 2Cu + a)Lod - aLot - {Cu + Lu){2Lot - 1) - {Li, Lq}, (3.38) 
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and two expressions for two components of the matrix Bq: 

B+ = KtBo+ = (1 - 2Lot)L2+ - (2Cit + a- l)Li+ + {1/t - 2L2t)Lo+, 



(3.39) 



Bt = ntBot = -Ll - {U, Lo} - 2L2tLot + L^t + L^t/t - Cu{Cu + a - 1) + do, (3.40) 

eq. (3.39) being another constraint, while eq. (3.40) just determines B^t in terms of Bqcl since 
it follows from eq. (3.28) that 

B = KtBo = B+a++Bda3+B^a- = (l-2Lot)L2-(2Cit+a-l)Li+(l/i-2L2t)Lo-do(73+aoa_. 

(3.41) 

Since B^t does not enter any other equations, eq. (3.40) has no other consequences. As for 
the Lax eqs. (3.23)-(3.25), for the following it is convenient to rewrite them in terms of just 
introduced matrix B, whose component B+ is determined by eq. (3.39) and the other two 
components Ba and S_ are still arbitrary: 

KtdtL2 = -Ktb_h2 + [L2, Ktb_hi - B], (3.42) 
KtdtU = Kt{b. - 6+)[L2, Lo] + [B, Li], (3.43) 

KtdtU = Ktb+Lo + [B - Ktb+Li, Lq]. (3.44) 

Thus, it is worth mentioning that, as a result of the above considerations, one is left with 4 
arbitrary functions of t: 6_, 6+, B^ and B_ entering the further analysis. The rest essentially 
depends on the way of resolving the constraints - eqs. (3.31)-(3.36). 

4 Different cases of constraint resolution 

Case I: Resolve eqs. (3.32), (3.34) as 

L2t = ^, Lot = I, (4.1) 
Then in place of eqs. (3.31) and (3.33) one gets, respectively, 

Ll^Ll, + Ll^^, (4.2) 

i^l-Ll, + Ll^\. (4.3) 

One then compares these with the consequences of anti-commutation of L2 with eq. (3.42) 
and Lo with eq. (3.44), respectively: 



13 



Ktdthl — —2K,tb^hl, 



(4.4) 



KtdtLl = 2Ktb+Ll, (4.5) 
from where it follows that two of the four arbitrary functions are not so anymore: 

tb^ = 1, b+ = 0. (4.6) 
Then the constraints (3.35) and (3.36) read: 

{k-2 + 2Cu + a)L2+ = 0, (4.7) 

(2Cit + a)Lo+ = 0. (4.8) 
Moreover, now the trace eqs. (3.37) and (3.38) become algebraic: 

(« - 2 + 2Cu + a)L2d = -{L2, Li} - (4-9) 



{2Cu + a)Lod = -{Li,Lo}-t (4-10) 



The eq. (3.41) now becomes 



B = -(2Cit + a - l)Li + 5^0-3 + S_(7_. (4.11) 
Case II: Resolve eqs. (3.32), (3.34) as 

L2t = |, ^0+ = 0, (4.12) 

considering Lot 7^ 1/2. Then eqs. (3.31) and (3.33) turn, respectively, into eq. (4.2) and 

{Lot + Lod){Lot + Lod - 1) = 0, (4.13) 

eq. (3.35) becomes 

(Ktb- -2 + 2Cu + a)L2+ = 0, (4.14) 

and eq. (3.36) gives 

Li+ = 0. (4.15) 

Taking also L2+ = would give a subcase of Case IV below, so we choose K,tb-—2+2Cit+a — 
in eq. (4.14). Eq. (3.37) becomes again algebraic, 

{L2,M = -^^|^, (4.16) 
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while eq. (3.38) remains intact so far. Eq. (3.39) simplifies to 

B+ = KtBo+ = (1 - 2Lot)L2+, (4.17) 

with right-hand side not identically zero. 

Then it is easy to verify that under the above conditions ODEs (3.43) and (3.44) can 
be solved only by taking Li = and Lq = 0. Now the relation k, — 2 + a = follows 
from eq. (4.16), and this is the equation (2.1) if one returns to the original a (recall that we 
redenoted a + 2//3 — a + l/K^a before). 

Eq. (4.2) and anticommutator of eq. (3.42) with L2 imply t6_ = 1 again and then Cu = 
follows from eqs. (4.14) and (4.16). Eq. (4.13) and Lq = give Lot(Lot — 1) = 0, and this 
together with other above facts entails Lot = as solution of eq. (3.38). Thus, only eq. (4.2) 
and eq. (3.42) together with eq. (4.17), 5+ = L2+, remain, containing two arbitrary functions, 
which ultimately leads to the Lax pair from Theorem 3, which makes sense for a < and any 
(nonnegative) value of /?. 

Case III: Resolve eqs. (3.32), (3.34) as 

L2+ = 0, Lot = \. (4.24) 

It turns out that Case III leads to a special type solution for which k = 2 and a = 0, which 
can be seen similarly to the above consideraton of Case II. This seems not very interesting 
since we have more general solution for /t = 2 and arbitrary admissible a. 
Case IV: Resolve eqs. (3.32), (3.34), (3.35) and (3.36) as 

L2+ = Li+ = Lo+ = 0, (4.25) 
Then eq. (3.39) gives B+ — 0, while eqs. (3.31) and (3.33) become, respectively, 

(L2t + L^d) {L2t + - l/t) = 0, (4.26) 

(Lot + Lod) (Lot + Lod - 1) = 0. (4.27) 

The '+ '-components of eqs. (3.42)-(3.44) become now trivial, while their diagonal components 
reduce to 

KtL'^^^ -Kth^{t)L2d, (4.28) 
Lid^Cu^ const., (4.29) 

KtL'^^ = Ktb+{t)Lod: (4.30) 

respectively. Adding trace ODE eq. (3.35) with eq. (4.28), taking into account eqs. (4.25) and 
(4.29), gives 
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Kt{L2t + L^dY = (2(Cu + Cirf) + a - 2)(L2t + L^a) - {Cu + Cu)/t, (4.31) 

and adding trace ODE eq. (3.36) with eq. (4.30), taking into account eqs. (4.25) and (4.29), 
gives 

Kt{Lot + LodY = -(2(Cii + Cid) + a){Lot + Lod) + (Cu + C^)- (4.32) 
However, eq. (4.27) implies that (Lq* + Lq^Y = 0, and eq. (4.32) becomes 

{2{Cu + Cm) + a){Lot + Lod) = Cu + C^- (4.33) 

Further consideration of the above constraints leads one to the conclusion that in this case 
there are only possibihties similar to Cases II and III above. 

5 Painleve equations for P = 2 and 4 

It is convenient to introduce 

^2 = tL2. (5.1) 



5.1 Case la 

Suppose 2Cit + a 7^ 0. Then eq. (4.8) gives Lo+ 
let 



= 0, therefore eq. (4.3) becomes Lq^ = i. So 



2-^od — Cod, Cod — ±1- (5.2) 

The "+"-component of eq. (3.44) now reads (2Cit + a - l)Li+ = 0. If Li+ = 0, then the 
"+" -component of eq. (3.43) gives L2+ — 0, and this reduces to a subcase of Case IV. So we 
take 

2Cu + a - 1 = 0. (5.3) 

Then eq. (4.7) becomes {k, — 1)L2+ = 0, and again taking L2+ = reduces to a subcase of 
Case IV since then the "+" -component of eq. (3.42) gives also Li+ = 0. Thus, we get k. = 1 
(i.e. l3 = 2). We introduce new variables: 

Li+ = 0, y^^-2^, l-^cPLo-, A_^cf>B_, A^^Ll^^(t>L,^. (5.4) 

The system of ODEs (3.42)-(3.44) is now 

tdtX2^[X2M-B], (5.5) 
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t%U^[X2M]+t[BMi (5 

tdtU^[B,U]. (5 

The diagonal component of eq. (5.7) is redundant, and its '+ '-component is trivial. The o 
nontrivial component of eq. (5.7) reads 

tdtl- = t^l- + CodA. - 2Bal- . (5 
The components of eqs. (5.5), (5.6) become 

tdtX^d = ^(A - - <f>X,_, (5. 

tdtX2+ = 2cf>X2d - 2{Lu - Bd)X2+, (5.5 

tdtX2- = 2(Lid - Bd)X2- - 2^(Ai - A_), (5.5 

t%L,a^^l_-tA_, (5. 

= -CodX2+ + 2i0Sd, (5.6 

i^^t^i = t^^Ai + Cod0X2_ - 2X2dZ_ - 2tiBaAi - L^A^). (5.6 
We now eliminate expressing it from eq. (5.6+) and using y from eq. (5.4), 

X2+ = -ca''y(t>- (5 

Substituting it into the constraint eqs. (4.2), (4.9) gives, respectively, 

4(^2, - Codt'y<pX2-) = 1, (5. 

2LidX2d - Codt'^yAi + 0^2- = (5. 
The remaining constraint to take into account, eq. (4.10), reduces to 

CodLid + 1-^ ~2~^' 

The ODE (5.6d) becomes, after substituting eq. (5.9) and using ^/-variable, 

tdtLid = -Codtyl- - A_. (5. 
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In fact the algebraic equations (5.10), (5.11) and (5.12) turn out to be first integrals of the 
remaining ODEs. Therefore it is possible and convenient to consider eqs. (5.5-) and (5.6-) as 
redundant. We next express I- from eq. (5.12) and substitute it and y-variable into eq. (5.8), 
which turns into 

tdtL,a = ty{Lu + (Co^a + l)/2) - A_, (5.14) 
which is identical to eq. (5.13) after ehminating /_ from the last by cq. (5.12). We express 
X2- from eq. (5.10), X2- = 40^2:^, plug it into eq. (5.11) and express Ai from there, 

Codt'yA, = ^^^^^ + 'iLuX^a + (5.15) 

After introducing 

r = -2CodX2d, Id = L^a + (Co^a + l)/2, (5.16) 
we substitute expressions for X2- and Ai into eq. (5.5d), and the last becomes 

- 1 

tr' = -2t''yA_ + — (2/^ - Co^a - l)r + Co<i(a - 1). (5.17) 

With the substitutions from eqs. (5.4), (5.9) and (5.16) done, eq. (5.5-I-) can be represented 
as 

tify)' = -iV + r - (2la - Co^a - l)t'y, (5.18) 

and eq. (5.14) now reads 

tl'^ = tyk-A^. (5.19) 

Thus, we are left with eqs. (5.17), (5.18) and (5.19). There is stil one arbitrary function, A_, 
left {Bd was already absorbed), but this can be effectively absorbed too. We differentiate 
eq. (5.18) to get 

t{t'y)" = -2t'yy' - 3tY + r' - 2t'yl', - {21, - Coaa){t'yy. (5.20) 

Then we subtract eq. (5.19) multiphed by 2t'^y from eq. (5.17), which cancels A^ there, and 
obtain 

t{r' - 2eyQ = -2iV + - (2/d - Co^a - l)r + Cod(a - 1). (5.21) 

Now we plug the right-hand sides of eq. (5.21) into eq. (5.20) multiplied by t: 

i2(i2y)''- -2tV-3iV-2iy+^^^-(2i<i-C'oda-l)r + C'od(a-l)-(2Zd-Coda)t(t2y)'. 

t y 

(5.22) 
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Next we put the expression for r from eq. (5.18) in eq. (5.22). Then all terms containing 
cancel and we obtain equation for y: 

II {y'f 1 ' , 3 C'odfl 2 , Cod(Q-i) 1 

This last equation is in fact a form of Painleve III. To bring it to the standard form of Painleve 
III, see e.g. [231 US], one just needs to define 



i = t-'/\ y = t-"'z = ez, 
then, in the new variables, eq. (5.23) is a standard Painleve III [231 ITU] : 

, ^ (£W _ + 4/ - + - '> - 1. (5.24) 

If one puts Cod = —1, then eq. (5.23) becomes exactly eq. (2.7a) from Theorem 1. 

5.2 Case lb 

Now take 2Cit + a = 0. Then eq. (4.8) is trivial. Using eq. (5.1), one obtains (k — 2)X2+ = 
from eq. (4.7). Eq. (3.41) now becomes B = Li^(r^ + 8^0-3 + B_a^, where functions Bd{t) and 
B_{t) are arbitrary. If /t 7^ 2, then X2+ = here. Therefore eq. (4.2) becomes X| = = j, 
while the '+'-component of eq. (3.42) reads simply (k — l)X2dLi+ = 0, which gives again 
K = 1, since taking Li+ = leads again to a subcase of Case IV. The further consideration is 
similar to the Case I-I. Now the system of ODEs (3.42)-(3.44) is 

tdtX2 = [X2M - B], (5.25) 

t''dtU = [X2M]+t[B,U], (5.26) 

tdtLo=[B,Lo]. (5.27) 

We again introduce new variables: 

Li+ = 0, y = ^-2 ^^'~/'"'\ /- = 0X2_, A_ = (t>B_, A = LL = 0i.i-. (5.28) 
(p t 

Let also 

2X2d = C2d, C2d = ±1, (5.29) 
then the only nontrivial '-'-component of eq. (5.25) reads 

U'_ = A. + 2{Lu - Bd)l- - C2d{Ai - A_). (5.30) 
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The components of eqs. (5.26), (5.27) become, respectively. 



t%Lu = -^l- + t{A, - A_), (5.26d) 
t''dt(l> = C2dLo+ + 2t(j){Bd - Lu) , (5.26+) 
t^dtAi = t^^A, + 2Lodl- - C2d(f>Lo- - 2t{BaAi - L^A^), (5.26-) 

tdtL^d = 0i^o- - (5.27d) 
tdtLQ+ = 2S<iLo+ - 20Lo<i, (5.27+) 

i9,Lo- = 2^A_ - 2SrfLo_, (5.27-) 

Three among six components of eqs. (5.26) and (5.27) are again redundant due to constraints- 
first integrals (4.3), (4.9) and (4.10). The algebraic equations (5.32), (5.33) and (5.34) below 
are first integrals of the remaining ODEs. Therefore it is possible and convenient to consider 
eqs. (5.26-) and (5.27-) as redundant. We ehminate Lo+, expressing it from eq. (5.26+) and 
using y from eq. (5.28), 

Lo+ = C2dt%(t). (5.31) 

Substituting it (and notations from eq. (5.28)) into the constraint eqs. (4.3), (4.10) gives, 
respectively, 

4(Lg, + C2di'#Lo-) = l, (5.32) 

2LuLm + (t>LQ- + C2dt'yA^ = -|. (5.33) 
The remaining constraint to consider, eq. (4.9) reads now 

C2dLu + l-^-"—^^^, (5.34) 
The ODE (5.26d) becomes, after substituting eq. (5.31) and using y-variable, 

tdtLu = -C2dtyl- + ^1 - (5.35) 

We next express /_ from eq. (5.34) and substitute it and variable y from eq. (5.4) into 
eq. (5.30), which turns into 

tdtLu = ty {lu + C2d ■ " ~ ^ ~ ^'" ^ + ^1 - (5.36) 
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which is identical to eq. (5.35) after ehminating l_ from the last by eq. (5.34). We express 



4L^ —1 

Lo- from eq. (5.32), Lq- = ~ 4C'."t^' (A ' then plug it into eq. (5.33) and express Ai from there, 



t'yA, = ^^1^ - 2C,,L,,Lo. - (5.37) 
Now we make the last change of variables, introduce 

r = -2C2dLo<i, Id = Lid + (C2<i(a - 1) - l)/2. (5.38) 

With the substitutions from eqs. (5.28), (5.31) and (5.38) performed, eq. (5.27+) transforms 
into 

t{t^yy = -tY + t^yi'lld - C2da + C2d + 1) + r, (5.39) 
When we substitute the expressions (5.31), (5.32) and (5.38) into eq. (5.27d), the last becomes 

tr' = ^^ + 2tV- (5.40) 
Next we plug expression (5.37) for Ai into eq. (5.36) and use (5.38) there to obtain 

t'yl'd - (tV + r)ld + ^ - + 1) + ^^^r - ^V- (5-41) 



Thus, we are left with eqs. (5.39), (5.40) and (5.41). There is still one arbitrary function, 
left {Bd was already absorbed), but this can be absorbed too. We differentiate eq. (5.39) to 
get 

t(fy)" = -2t^yy' - 3tY + r' + 2t^yl'^ + (t'i/)'(2/rf - Ca.a + C2d). (5.42) 

Then we add eq. (5.40) and eq. (5.41) multiplied by 2, which cancels terms with A_ there, 
and obtain 

t{r' + 2fyQ = + 2{ey^ + r)ld - C^da^r + 1) + {C^d + l)r. (5.43) 

Now we plug the right-hand sides of eq. (5.43) into eq. (5.42) multiphed by t: 

t\ey)'' = -2t''yy'-2,ey^+''^^+2{tY+r)ld-C2da{r+l)^^^ 

(5.44) 

We then express Id from eq. (5.39) and put it into eq. (5.44). All the terms containing r then 
happily cancel and an ODE for y results: 

/=(|^-i,' + ,3 + ^i(|^,._^_J_. (5.45) 
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This last equation is in fact a form of Painleve III. To bring it to the standard form of Painleve 
III, see e.g. US], one just needs to define 



e = t-''\ y = t-^'^z = 
then, in the new variables, eq. (5.45) is a standard Painleve III [23| 119]: 

= (fW - + 4.^ + _ 1 _ 1^. ,5,46) 

Z ^ ^ Z ^ 

If one puts C2d = — 1, then eq. (5.45) becomes exactly eq. (2.7b) from Theorem 1. 
5.3 Case Ic 

Also take 2Cit + a = but resolve eq. (3.35) as k = 2, then eqs. (3.37) and (3.38) become, 
respectively, 

{X,M} = ~, (5.47) 

{Li,Lo} = -|. (5.48) 

Take now the two arbitray functions in eq. (3.41) to be do = = 0. Then eq. (3.41) 
simplifies to B = Li. The triple of ODEs (3.42)-(3.44) now reads, when also taking into 
account eq. (5.1), 

2tdtX2 = [X2,Li], t^dtU = [X2MI 2tatLo = [Li,Lo]. (5.49) 
Let us now change variables since we want to use some symmetry between X2 and Lg, i.e. let 

S = X2 + Lo, A = X2- Lq. (5.50) 

One can rewrite all the remaining equations in terms of Li, S, and A. Recall also constraint 
eqs. (4.2) and (4.3). Their sum and difference give, respectively, 

S^ + A^ = l, {S,A} = 0, (5.51) 
while the sum and difference of eqs. (5.47) and (5.48) give 

{S, Li} = -a, {Li, A} = 0. (5.52) 
At last, the system of ODEs (5.49) is equivalent to the following system: 

2tdtS = [A,U], (5.53) 
2t''dtU = -[S,Al (5.54) 



22 



2tdtA^[S,U]. (5.55) 
Next we rewrite the system in diagonal-anti-diagonal form (S — S^cts + Sa etc.): 

Sj + Si + Al + Al^l, 2SdAd + {Sa, Aa} = 0, (5.51) 

2SdL,d + {Sa, L,a} = -a, 2LuAd + ^ J = 0, (5.52) 

2tdtSa=[Aa,L,a]az, {5.53d) 



tdtSa = Ada^Lia - Lua^Aa, (5.53a) 

2edtLu = -[Sa.Aa]a^, {5.5Ad) 

t^dtLia = A^a^Sa - Sda^Aa, (5.54a) 

2tdtAd^[Sa,L,a], (5.55rf) 

tdtAa = SdCTsLia - L^a^Sa (5.55a) 

It turns out that likely a simplest way to find the solution here is to put Aa — 0. Then the 
system is greatly simplified and its non-trivial solution can be found by choosing also 

Sd = Lid = 0. 

Then the second equations in eqs. (5.51), (5.52), as well as eqs. (5.53d), (5.54d) and (5.55a) 
are trivially satisfied and what remains is 

Si + Al = l, (5.56) 

{Sa,Lia}^-a, (5.57) 

tdtSa = AdasLia, (5.58) 

t'^dtLia = Ada^Sa, (5.59) 

2tdtAd^[Sa,Lia]. (5.60) 

Now, after expressing Sa from eq. (5.59) and plugging it into eqs. (5.56), (5.57) and (5.60), 
they turn into, respectively. 
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t\d,L^^y = AliAl-l), 



(5.61) 



t^[Lia,dtLia](r3 = aAd, (5.62) 

d,iAl)=td,iLl). (5.63) 

Differentiating e.g. eq. (5.57) and using eq. (5.58) one can see ttiat eq. (5.58) can be considered 
as redundant - it is a consequence of tfie otliers. To close tlie system of tfie last tliree written 
out equations, it is convenient to use tlie identity: 

{Li,, dtL,,y - dtL.a]' = ALlidtLia)'. (5.64) 
We substitute eqs. (5.61), (5.63) and the square of eq. (5.62) into eq. (5.64) and obtain: 

A{Al - l)Ll, = AitdtA^y - a\ (5.65) 
Next we differentiate eq. (5.65) and find: 

A{Al - l)dtL\^ + ^A^dtAA = KtdtAd)dt{tdtAd). (5.66) 

Finally we substitute from eq. (5.65) and dtL\^ from eq. (5.63) into eq. (5.66), which 
yields an ODE for the single variable Ad. 

t{Al - l){tA,)' = A,{tA',r + ^-A',{Al - 2) + Q - ^) A,, (5.67) 

where again ' = dt- One can now recognize the last equation as almost identical to the equation 
(1.16) of the paper [32] about the hard edge of unitary- invariant random matrix ensembles 
(there /3 = 2) and Bessel kernel. By changing the independent variable in eq. (5.67) to 
^ = 4/t, the last becomes exactly that equation from [32], which is shown there to be a simple 
transformation of a sigma- form of Painleve III: setting A^/A = d{^R{^))/d^, one can show 
that -R(0 satisfies a sigma- form of Painleve III [32], i.e. this ^R{^) is a particular classical 
Painleve III Hamiltonian. 

Thus, for /3 = 4 hard edge we also have a Painleve III representation. 

6 Explicit Lax Pairs 
6.1 Case la 

Here (3 = 2, Lo+ = 0, Lg^ = 1/4 and Cu = (1 — a)/2. Anticipating a simpler result, we put 
Bfi = Ba = and Lq^ = —1/2. From the '-'-component of eq. (5.7) it follows that Lq- = const 
and we put Lq- = 0, in order to get the simplest solution. Then we have X2d = r/2 from 
eq. (5.16), and eq. (5.12) gives 
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^" 2 ' ^ ' 

Clearly now eq. (5.6d) is satisfied by eq. (5.12). It is convenient to express the antidiagonal 
components from eq. (5.6), = ^^crs-^'ia- Substituting the last expression into eqs. (4.2), 
(4.9) and (5.5d) yields, respectively, 

t'{L[J' = lir' - 1), (6.2) 
t'[L^^,L'Ja,^^{r + l), (6.3) 

2t{Liy = r'. (6.4) 
Now we use the identity eq. (5.64) and plug the last three equations into it to derive 

4(r2 - l)Ll^ = {tr'f - (a - lf{r + 1)^. (6.5) 

This is to be used together with eq. (6.4) to eliminate and obtain an ODE for r. We 
differentiate eq. (6.5) and, comparing with eq. (6.4), we obtain the equation satisfied by r, 

t(r' - l){tr')' = r{tr'f + -(r' - if - (a - lf{r + if. (6.6) 
If we introduce function q such that 

r = -l + 2g^ (6.7) 
which is inspired by the hard-to-soft edge limit, see below, then the new q will satisfy 

t{q^ - l){tq'y = q{tq'f + \{q^ - if - ^^^g, (6.8) 

which is nothing but an instance of the Tracy- Widom equation (5.67) (up to change a — I a) 
for /3 = 2 hard edge. 

Now we write out the explicit Lax pair for this case using the last findings and the previous 
consideration ending in a standard Painleve III for y — L'i_,_/Li+, see eq. (5.23) (one should 
put Cod = — 1 there to be consistent with the current consideration). First we recall that 
X2 = and so L2d = X2d/t = r/{2t) and L2+ = ty(j) follows from eq. (5.9). Next, from 
eqs. (5. 9), (5. 4) and (5.12), an explicit connection between current r and y is 

r^t^(y' + y'') + (a + l)t^y. (6.9) 
Recalling eq. (6.5) and the fact that = Li_,_Li_ = 4>Li_, we find the expression for Li_: 

^'- = 5^- '»-^-2(<-i)(^ + i)- ("0) 
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Prom eq. (5.10) we find 



- ^'-1 

Finally the Lax matrices read: 



(6.11) 



\ r f A r f ^ ^4J(I) 



n 1 / ' - I h 1_Q I) ^2 - \ 1-r 

and, using eqs. (3.21), (3.22), S = and (3.40), we get 



Thus, the whole Lax pair reads: 



l+r ^ _|_ 

2tx2 „ a; a;2 

1 1 1-r I 1-a 



lul. _|_ in" _|_ 1 



where y = (f)'/(/) is a solution of Painleve 111 eq. (2.7a), h is defined in eq. (6.10) and r is a 
solution of eq. (6.6), which is related with y by eq. (6.9). Also, y can be expressed in terms 
of r from eq. (5.17) (where ld — now) as 

^2 _ I 

^ " t2(tr' + (a-l)(r + l))" ^^'"^^^ 

With r replaced by g according to eq. (6.7), the above Lax pair will take exactly the form of 
the first Lax pair from Theorem 1. 

6.2 Case lb 

Also (5 = 2 here. Since L2+ = 0, eq. (4.2) gives = 1/(4^^). To get a simpler solution, we 
choose — —1, i-e. — — ^- Prom eq. (5.31) it follows then that Lo+ = — = —t^y(t>- 
The simplest solution is obtained if one takes S = Li in eqs. (5.25)-(5.27). Then from 
eq. (5.28) we find that (j)' = ycf) and = Ai. It is also the simplest to take Id = —C2dl- = 0, 
which solves eq. (5.36) and means that L2- = 0, by eqs. (5.28) and (5.1). From eq. (5.34) one 
finds then that Lid — |- From eq. (5.38) it follows now that Lq^ = r/2, from eq. (5.32) one 
finds 

- 1 
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while eq. (5.37), together with eqs. (5.34) and (5.28), leads to 

^'- = T = 5^' "-^ + ^. + 1). (6.15) 

Now we recall that Bi = by eq. (3.22), B_i = -dtL2t ■ I + b-{t)L2 = 1/(2^^) . / + La/t by 
eq. (3.21). Finally, from eq. (3.40), recalling also that here — Cu — —a/2, we find 

1 + r h a 
^ ^ " 4^ ~ 2- 
Gathering all this information together, we can write down 

00\ r _f 0\ r _f '-y" -i'y<t^ 

n 1 J ' -^1 — I h _n J ' -^0—1 r^-l 1-r 
^ t J \ 4:t^y<P " / \ 4*22/,/. 2 



B -B A+^I _ A + r h ^\ r,U I. ^ 

Bo-Bo^-I + -U-^-^-^--j.I+- 



2 

h a 



and the whole Lax pair reads: 



h 1 _|_ r^ — 1 _L ^ -I- 

4t2j/</> X 4t2y<)!) tx2 X 2 



V 2^2 4t3yy \^ ^ W^-J J 

where y — (f)' / (f) \s di, solution of Painleve III eq. (5.45) (with C2d — —1), h is given in eq. (6.15), 
y and r arc tied by the system, see eqs. (5.39) and (5.40), where now la — ^ and A_ — Ai is 
given by eq. (6.15), 

t(t'2/)' = r-tV + atV (6.17) 
- 1 

tr' ^—— + a{r + l). (6.18) 
t^y 

It is easy to find from the above that r itself satisfies the following second-order ODE: 

t{r^ - l){tr')' = r{tr'f + -(r^ - 1)^ - a'ir + 1)1 (6.19) 

If we again introduce function q as in eq. (6.7), r — —1 -\- 2q^, which is inspired by the 
hard-to-soft edge limit, see below, then the new q will satisfy 

t{q'-mq')'^q{tq'f + \{q'-lf-^^q, 
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which is nothing but the Tracy- Widom equation (2.12) (or (5.67)). Then the above Lax pair 
will take exactly the form of the second Lax pair from Theorem 1. 



6.3 Case Ic 

Here (3 — A, one has Li^ — 0, Sa — tL2d + Lq^ — and — tL2a — -^oa = by construction. 
Let Ad = q - the solution of hard edge Tracy-Widom eq. (5.67), then immediately from 
eqs. (5.50), (5.1) and the fact that Sd — we get 

L2d = |, Lod = -|. (6.20) 

One has also Sa = tL2a + L^a = '2fL2„ = 2Loa- Eq- (5.56) yields = S+S_ = 1 — q"^, while 
eq. (5.65) gives Lf^ = Li_|_Li_ = "^4(g2_i) • From eqs. (5.57) and (5.60) one finds 

2S+Li^ = 2tq' - a, 2S-Li+ = -2tq' - a, (6.21) 
which, together with eqs. (5.56) and (5.65) gives also the ratios 

Li+ 2tq' + a Li_ 2tq' - a 



S+ 2(g2-i)' S- 2(g2_i)- 
We recall also the eqs. (5.58) and (5.59) which mean, respectively, that 



(6.22) 



tS'_^ = qLi+, tS'_ = (6.23) 

t'^L[^ = qS+, t^L\_ = -qS-. (6.24) 
Prom eqs. (6.22) and (6.23) we express the logarithmic derivatives of and -S"-, 

S'+ ^ qq' aq ^ ^ /^^ 

S+ q^-l 2t{q^-iy S_ q^-l 2t{q^ - 1)' ^ ' 

and the components Li+ and Li_ are then expressed from eq. (6.22), which completes finding 
of the components of matrix L in terms of the solution q of eq. (5.67). Recalling also that 
Lit — 1/(2^), Lit — —a/2 and Lot — 1/2, we can write out 



/ 1+9 








l-q I ) 




V 2t 


2t / 





a 



2tq'+a ^ \ / i_q S+_ 



2 2(g2-l)"^+ \ T — i ~Y' ~ 

2tq'-a c _a I' ^0 - I ^- 1+q j ) 

"2(g2-l)'-'- 2 / V 2 2 

where the functions S+ and S_ are completely determined by eqs. (2.13). Now we recall that 
5i = by eq. (3.22), 5_i = 1/(21'^) ■ I + Li/t by eq. (3.21), Bo = U/{2t) by eqs. (3.41) and 
B = U. Finally, by eq. (3.40), 

{2tq'r - a' q\ a{a - 2) 

= - 4(,^ - 1) + T + 

Thus, after gathering all this together, the Lax pair for ^ = 4 is the one of Theorem 2. 
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6.4 Case II 



This is a case where k, = 2 — a, induced by requiring that L2t = 1/(2^) and Lot 7^ V^- The 
final result depends on two arbitrary functions, which we chose in such a way that the Lax pair 
bears some resemblance with the previous Lax pairs where all the functions are determined, 
unlike here. 

Again let X2 = tL2 as in eq. (5.1). One of the components of eq. (3.42) is redundant 
due to eq. (4.2) being its first integral. The other two components depend on two arbitrary 
functions as e.g. 

i^t^2d — ~b-^2+, KtX2^ — 2bclX2+, 
where, using eqs. (4.17) and (4.23) we represented the matrix B from eq. (3.41) as 

We denote the two arbitray components as — r/2 and X2+ — (j). Then eq. (4.2) gives 
X2- — —^-4^, and the other components of matrix L are already known from section 4. Using 
eqs. (3.21), (3.22) and (3.40) we find the non-arbitrary elements of matrix B. The final result 
is the Lax pair of Theorem 3. 

One can easily verify that, for any two functions r{t) and 0(t), the zero curvature equations 
are satisfied identically with this pair. The Fokker-Planck equation is satisfied by the first 
component of the eigenvector of the pair provided that k — 2 — a only, but with no other 
restrictions. However, the above Lax pair means that 

^ 2^ + + ^'^^^ 

and 

which are multiples of each other, and so this J-"^ in fact also satisfies a much simpler equation, 

KtdtT^ — x{x + K)dxT^ , 
and, as a consequence of this and the FP eq. (2.2), also satisfies 

both of which can be solved and with proper initial conditions for J^^ result in eqs. (2.16) 
and (2.17) of Theorem 3 for J^" and G. 
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7 Hard edge to soft edge limit transition 



This transition was considered before in [TO] in terms of Fredholm determinants and its 
kernels and in terms of stochastic 1st order ODEs in ^] . Here we track how it goes in terms 
of the derived above Painleve equations, which partly complements considerations of [18], and 
related Lax pairs, for the cases /3 = 2 and 4. The transition occurs in the limit as a — )■ oo, 
then, if we denote a = f , the independent variables x and t in the FP equation (2.2) scale as 
m eq. (2.18), i.e. x^a{l + a-^^^x), t ^-^(l + a-'^^H). Then, since the partial derivatives 
scale as 

O^^a-'/'d,, dt-^a'/%, (7.1) 
the hard edge FP operator Vh scales as 

Vh = (Ktdt + x^d,, + (ax -x^- l/t)d,) ^ a^^^Vs = a^^^Kdt + d,, + (t - x^)d^), 

where we have the soft edge FP operator written out by [S] on the right-hand side. This simply 
and clearly justifies the name of the transition, and shows also that the solution of eq. (2.2), 
J-'^ , scales as J-'^ — )■ a~'^^^J-'^ , with J'^ being the solution of soft edge FP equation [8]. 

7.1 The transition for (3 = 2 

We only write down the derivation for the Case la of Lax pair (2.6a) since the transition for 
Case lib of pair (2.6b) gives the same final result (even though the scaling of functios involved 
is different which can be most easily figured out from the algebraic identities presented in the 
remark after Theorem 1, e.g. the scaling of y changes to y a^{l + a~^^^y) instead of that in 
eq. (7.15) below and G corresponding to eq. (7.19) below changes to G — >• e" however, 
the scaling of q does not change). Recall the defining system of first-order ODE eqs. (6.9) and 
(6.13) obtained in section 6.1 for functions denoted y and r {y being the solution of Painleve 
III eq. (2.7a)): 

t'y' = r - t^y^ - {a + l)t'y, (7.2) 

t^yr' = (r + 1) (r - 1 - (a - l)t^y) . (7.3) 
We consider the scaling limit of these equations, using eqs. (2.18) and (7.1) and substituting 

y yo{a) + a^^i/i, r ro(a) + a^Ti. (7.4) 
Then the highest order in a gives, repectively, 

ro - 2a-^yo - aSl = 0, (ro + 1) (ro - 1 - 2a-^y^) = 0, 

with possible solutions ro = — 1, yo = — «^ and tq = 3, yo = a^. Taking the next terms in a 
shows that only the first solution is consistent. Then the next order terms in (7.2) give 
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Similarly, the next order terms in eq. (7.3) yield 

= ri(2a^^-^2/i - a^Ti). (7.6) 
Equations (7.5) and (7.6) consistently imply that the scaling dimensions are 

A, = 1, Ar = -| (7.7) 
and thus, eqs. (7.5) and (7.6) finally become the corresponding system for the soft edge limit: 

y[ = r,-yl + t, (7.8) 

r[ = 2yin. (7.9) 
Solving it for yi and ri, respectively, one gets 

y'l = 2yl-2ty,-l (7.10) 

and 

2rir'l = {r[y + Atrl + 4r^ (7.11) 

Equations (7.10) and (7.11) can, of course, also be obtained by taking the scaling hmit of the 
corresponding second-order ODEs eq. (2.7a) and eq. (6.6). The transformation 

yi = I'l'y. t = 

brings eq. (7.12) to the standard Painleve II form: 

y"{s) = sy + 2y'-^, (7.12) 
while (7.11) is a form of Painleve XXXIV equation, see e.g. [12]. If one substitutes 

ri = 2q^, (7.13) 
then one obtains the Painleve II equation with free parameter zero: 

g"(t) = tg + 2g3, (7.14) 

which has the Hastings-McLeod solution relevant to the soft edge limit pJLj. It may seem 
surprising that another Painleve II, with free parameter —1/2 also arises here. This, however, 
has a natural explanation in the fact that the solutions of the two Painleve II equations are 
tied by the Gambler relation, see e.g. recent nice Painleve II review 



31 



where means a solution of Painleve II in the standard form with free parameter i/, and 
e = ±1. This relation also naturally explains the discrepancy of independent variables in 
eqs. (7.12) and (7.14). 

Now let us consider the soft edge limit for the hard edge Lax pair (6.12) of subsection 6.1 
(this is the Lax pair (2.6a) in terms of r instead of q from cq. (2.12)). 

One should have L — )■ ct'^I^Ls and B — )■ a^^^Bs, according to eq. (7.2), where Ls and 
Bs are the limiting soft edge matrices to be found. One finds, however, with the help of 
eqs. (2.18), (7.1) and 



y ^ -a^{l - a~'^/^y), r ^- -1 + q;"^/^, (7.15) 



that 



^2/3„ai/3t _ r{x-y) ^2 _ f _ r 



(7.16) 



the exponential factors appear because scales as 



^ e-"'^'V, (7.18) 

and new again satisfies the relation 0' = y(j) with the new y from (7.15). Rescaling J-"^ 
according to (2.22) removes the additional a^^^^ factors in the non-diagonal entries of L. The 
exponential factors e^"'^^^* as well as the term a^/^ in the last entry of B in (7.17) cancel upon 
rescaling of the second component of the eigenvector G as 

G ^ e"'^'*^^, (7.19) 

which determines the correct scaling of G in the transition to the soft edge. Thus, finally we 
obtain the corresponding soft edge limit Lax pair: 

d.( ^l)=Ls( ^l), dt( ^l)=Bs( (7-20) 

where 



I - ^5 I , G^ ~^'\ G^ 
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and q is the Hastings-McLeod solution of the Painleve II equation (7.14) and r = 2q^ for the 
rescaled r on the right-hand side of the second relation in (7.15), see eq. (7.13). Now we infer 
from eqs. (7.9) and (7.13) that our turns into q and we obtain the Lax pair of Theorem 4. 



7.2 The transition for /3 = 4 

Here it is convenient to consider the scaling limit of the hard edge Tracy- Widom equation 
(2.12). It turns out that under the scaling 

q -> a'^'^q (7.23) 

and the scalings from (2.18) and (7.1) for t and dt this equation turns exactly into the Painleve 
II eq. (7.14) with free parameter zero, q" = tq + 2q^. Thus, eq. (7.23) gives the right scaling 
for g, and then equations (2.13) determine that functions 5"+ and scale as ~ 1 in a. Thus, 
in the limit we indeed get eq. (2.24). 

Consider now the limit for the Lax pair involving the above functions. This time L rescales 
correctly from the beginning: 



(7.24) 




x—u—q S+x 

F ~xT+, ), (7.25) 
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where u is the function satisfying eqs. (2.21) in terms of the new q satisfying eq. (7.14). 
Therefore we infer also that here both of the eigenvector components rescale as in eq. (2.25). 



8 Conclusions 

We established the correspondence of solutions to FP equation for hard-edge probabilities of 
spiked large /3-Wishart (/3-Laguerre) matrices for special values of /3 = 2, 4 with eigenvectors of 
certain Lax pairs for Painleve III. Also it is shown that there is no such simple correspondence 
for other values of /3, as is the case with the FP equation for the corresponding soft-edge 
probabilities [S]. The only difference is the case when a special relation eq. (2.1) holds between 
the additional parameter a here and (3. Then the bounded solution of the hard-edge FP 
equation is explicit elementary function, eq. (2.16), which can be seen as a one-parameter 
generalization of the Gumbel distribution for extreme values. 
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Thus, the relation of the FP solutions, which represent non-stationary euclidean time 
eigenfunctions for quantum Painleve Hamiltonians, with /3 playing the role of the Planck con- 
stant (or its inverse), see |25], with the classical Painleve equations and classical integrability 
in general remains an interesting open problem, further considered in [30] . 
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